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ABSTRACT 



We study the effects of intrinsic transverse momentum smearing on J/ip production 
both at the Tevatron and at HERA. For the case of large-pr J /i^ production at the 
Tevatron, the effects due to smearing are mild. On the other hand, inelastic J /if) 
photoproduction at HERA is very sensitive to the kr smearing and, in fact, with a 
reasonable value of {kx) it is possible to resolve the large- 2; discrepancy seen by com- 
paring non-relativistic QCD (NRQCD) predictions with the HERA data. We conclude 
that, with the present kinematic cuts, photoproduction at HERA is not a good test of 
NRQCD. 
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Over the last few years, there has been a considerable advance in the understanding 
of quarkonium physics due to the development of the non-relativistic effective field 
theory of QCD, called non-relativistic QCD (NRQCD) The Lagrangian for this 
effective theory is obtained from the full QCD Lagrangian by neglecting all states with 
momenta larger than a cutoff of the order of the heavy quark mass, m, and accounting 
for this exclusion by introducing new interactions in the effective Lagrangian, which 
are local since the excluded states are relativistic. Beyond the leading order in 1/m 
the effective theory is non-renormalisable. The scale m is an ultraviolet cut-off for the 
physics of the bound state; however the latter is more intimately tied to the scales 
mv and mv^, where v is the relative velocity of the quarks in the bound state. The 
physical quarkonium state admits of a Fock expansion in v, and it turns out that the 
QQ states appear in either colour- singlet or colour-octet configurations in this series. 
Of course the physical state must be a colour- singlet, so that a colour-octet QQ state 
is connected to the physical state by the emission of one or more soft gluons. In 
spite of the non-perturbative nature of the soft gluon emissions, the effective theory 
still gives useful information about the intermediate octet states. This is because the 
dominant transitions that occur from colour-octet to physical colour-singlet states are 
via El or Ml transitions with higher multipoles being suppressed by powers of v. It 
then becomes possible to use the usual selection rules for these radiative transitions 
to keep account of the quantum numbers of the octet states, so that the production 
of a QQ pair in a octet state can be calculated and its transition to a physical singlet 
state can be specified by a non-perturbative matrix element. The cross-section for the 
production of a meson H then takes on the following factorised form: 

= E ^.(oSC'^'Lj)), (1) 

n={a,S,L,J} '"' 

where the F^s are the short-distance coefficients and the On are local 4-fermion oper- 
ators, of naive dimension dn, describing the long-distance physics. The short- distance 
coefficients are associated with the production of a QQ pair with the colour and angu- 
lar momentum quantum numbers indexed by n. These involve momenta of the order 
of m or larger and can be calculated in a perturbation expansion in the QCD cou- 
pling as{m). The QQ pair so produced has a separation of the order of 1/m which 
is pointlike on the scale of the quarkonium wavefunction, which is of order l/{mv). 
The non-perturbative long-distance factor (O^) is proportional to the probability for 
a pointlike QQ pair in the state n to form a bound state H. 

The existence of the colour-octet components of the quarkonium wave function is 
the new feature of the NRQCD approach. Before the development of NRQCD, the 
production and decay of quarkonia were treated within the framework of the colour- 
singlet model @, 1^. In this model, it is assumed that the QQ pair is formed in the 
short-distance process in a colour-singlet state. The corrections from terms higher 
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order in v were neglected. While this model gave a reasonable description of low- 
energy J/ip data, it was known that it was incomplete because of an inconsistency in 
the treatment of the P-state quarkonia. This was due to a non-factorising infra-red 
divergence, noted first in the application of the colour-singlet model to Xc decays 0, 
and the proper resolution of this problem was obtained only by including the colour- 
octet components in the treatment of the P-states 0. The colour-octet components, 
however, had a more dramatic impact on the phenomenology of P-state charmonium 
production at large px at the Tevatron pp collider |^ where the colour-singlet model 
was seen to fail miserably. While the inclusion of the colour-octet components for the 
P-states was necessary from the requirement of theoretical consistency, there was no 
such problem with the S states because the corresponding amplitude was finite and the 
colour-octet components were suppressed compared to the colour-singlet component 
by 0{v'^). But the data on direct J/ip and ip' production at the Tevatron [0 seem to 
indicate an important contribution from the colour-octet components for the S'-states 
as well [§. 

While it is clear that the correct description of the Tevatron large-pr data requires 
that the colour-octet components of the quarkonium wave function have to be taken 
into account, the major problem is that the corresponding long-distance matrix ele- 
ments are a priori unknown and can be obtained only by fitting to the Tevatron data 
Q. Clearly it is important to have other tests of NRQCD, so that the matrix elements 
which are obtained from the Tevatron data can be determined independently in other 
experiments. 

One important cross-check is the inelastic photoproduction of J/ if) at the HERA ep 



collider |jTO[. The inelasticity of the events is ensured by choosing z = pp ■ pj/^/pp ■ p^ 
to be sufficiently smaller than one. An upper limit of z ~ 0.8 would seem to be 
an appropriate choice An additional cut needs to be imposed on the pt of the 



J/ip to ensure that the production process occurs at short distance. In the HERA 
experiments, the z distributions have been studied using px > I GeV. The surprising 
feature of the comparisons |l3l of the NRQCD results with the data from HERA is 
that the colour-singlet model prediction is in agreement with the data while including 
the colour-octet component leads to violent disagreement with the data at large z. 
While the colour-singlet cross section dominates in most of the low-z region, the colour- 
octet contribution increases steeply in the large- 2; (0.8 < z < 0.9) region and this rise 
is not seen in the data. In these comparisons, the values of the non-perturbative 
matrix elements are taken to be those determined from a fit to the Tevatron large- 
ly data. Naively, one would think that this points to a failure of NRQCD. But this 
conclusion is premature and, in our opinion, incorrect. The reason is that while at the 
Tevatron the measured pt of the J/ip is greater than about 5 GeV, at HERA the pt 
can be as small as 0{1) GeV. At such small values of pt (and also for z very close 
to unity), there could be significant perturbative and non-perturbative soft physics 
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effects. One way to parametrise tliese effects is to include tlie effects of tlie transverse 
momentum smearing of tlie partons inside the proton. This is the purpose of the 
present paper. We study the effects of the parton transverse momentum, /ct, on the 
J /if) distributions both at the Tevatron and at HERA. We demonstrate that the z 
distribution measured at HERA is particularly sensitive to the effects of kx smearing, 
and argue that inelastic photoproduction at HERA, with the present kinematic cuts, 
is not a clean test of NRQCD. Other effects such as soft-gluon resummation and 



the breakdown of NRQCD factorisation near z = 1 [|T5| have been discussed in the 
context of this discrepancy. 

In our work, we first consider the effects of kx smearing on the extraction of the 
non-perturbative matrix elements from the large-pr direct J/ip production data from 
the CDF experiment at the Tevatron. Before we discuss the fits, we recall the inputs 
that go into the theoretical computations of Ref. where these fits are performed in 
the absence of smearing. The direct J/ip production (i.e. S'-state production, with the 
P-wave contributions removed) cross section in the NRQCD approach receives con- 
tributions from 1) the colour- singlet '^S'^^ channel, 2) the colour-octet ^Pf^ connected 
to the physical J/ip state by an El transition, 3) the colour-octet ^Sq'^ channel which 
makes a Ml transition to the physical state, and 4) the colour-octet '^S'f channel 
which connects to the J/ip state by a double El transition. All three colour-octet 
channels contribute at O(f^). The non-perturbative parameter for the colour- singlet 
channel (i.e. the radial wave-function at the origin) can be extracted from J/ip lep- 
tonic decay or estimated from potential model calculations. Given this input, the three 
non-perturbative parameters (O('^Pf^)), {0{^Sq^)), (O('^S'f')) (which we call matrix 
elements Mi, M2 and M3 respectively) have to be extracted by fitting to the CDF 
data. It turns out that for px values greater than about 4 GeV, the pt dependence 
of the short-distance coefficients corresponding to the ^-Pj^' and the ^Sq channels are 
identical. The ^S'f channel on the other hand has a different px distribution: the rea- 
son for this is that a fragmentation-type diagram where a single gluon attaches itself 
to the cc pair is present only for this channel and not for the ^Pj and the ^S^^ chan- 
nels. Due to the existence of this fragmentation-type diagram, the contribution from 
the ^S'f channel dominates at large pt and so it is possible to use the experimental 
data points in this pt region to fit the value of the corresponding non-perturbative 
matrix element M3. On the other hand, the contributions due to the ^Pp and the 
^^o*^^ channels are dominant at low p^ and one can use data the points at the lower 
end of the measured pt spectrum to fit a linear combination of Mi and M2. The linear 
combination turns out to be Mi/ ml + M2/3. We use the phrase 'fragmentation-type' 
guardedly. The subprocess cross sections to which we refer here are obtained from 
fixed-order perturbation theory calculations. In a genuine fragmentation calculation, 
the effects of Altarelli-Parisi evolution of the fragmentation functions are taken into 
account, and these go beyond fixed-order perturbation theory. In comparing the ^S'f^ 
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Figure 1: The CDF data for Bda/dpT (in nb/GeV) for J/ip production at 1.8 TeV 
with —0.6 < T] < 0.6, compared to the model predictions with (kr) = 0, 0.7, 1.0 GeV, 
respectively. 



prediction to the large-p-r data, one must correct for the effect of this evolution, and 
the prescription for this correction (as given in Ref. 0]) is to muhiply the number 
obtained from the fixed order calculation with the correction factor R, where R is the 
ratio of the cross sections (computed using the fragmentation approach) at the scales 
Pt and rric i.e. with and without the evolution of the fragmentation functions. 

For the Tevatron, we perform similar fits to those described above except that 
we now include the effects of intrinsic transverse momentum smearing. The inclusion 
of kx smearing in large-pr reactions is done in the standard fashion, see for example 
Ref. Each incoming parton is given an intrinsic transverse momentum, with a 



gaussian distribution whose width is an adjustable parameter. The gaussian shape is 
motivated by measurements of the transverse momentum distribution of , 77 and 
TTTT pairs produced in hadron-hadron collisions. The cross section with the inclusion of 



4 



smearing is therefore given by 

E^j^ipp^ J/i^X) = ^ j (f'kxa j (f'kTb j dxi J dx2 



a,b,c=g,g ' 



1 d& 

/a/p(Xl, kra, Q^)fb/pix2, kTb, Q^)--j{ab ^ 25+1^^^)^^^ + £ + u - M2), (2) 

IT dt 



where 



^2^ _ Fa/p{Xi,Q^)D{kTa) 



and similarly for fb/p{x2, kxb, Q^)- Fa/p is the usual parton distribution, and the func- 
tion Dikx) is the intrinsic transverse momentum distribution 

D{kT) = -exp(-624), (4) 

TT 

with 



and normalised such that 



(kr) = ^, (5) 



j d^krOikr) = 1- (6) 
For the numerical phase-space integration, it is convenient to use the change of variables 



1^^ 



kxi = K ( In— j , {i = a, h) (7) 



so that 

.2 



/ kndkxi dqi — . (8) 

JO Jo 2qi 

In terms of the variable Qi the Gaussian distribution becomes 

Dig.) = - {qf''" . (9) 
vr 

By choosing k such that is slightly larger than unity, we obtain a relatively smooth 
integrand. 

The results of the effects of including intrinsic kx for J / ip production at the Teva- 
tron are shown in Fig. 1, where we consider the J/ip px distribution, Bda/dpr, with 
B the branching ratio of the J/ip into ^~^fi~. The cuts on the pseudorapidity of the 
J/ip are those used by the CDF experiment i.e. —0.6 < rj < 0.6. In our computations 
we use MRS(D'_) parton densities |0 with factorisation scale Q = Mt = J Ml + p"^. 
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We present the results for three different values of (/ct), viz. {kr) = 0, 0.7, 1.0 GeV. 
The choice (fc^) = corresponds to no smearing. For the colour-singlet parameter, 
(C(^5'f' )), we use the value 0.81 GeV^. This corresponds, of course, to the case where 
there is no smearing. For the case when smearing is included we need to know the effect 
of the smearing on the value of this parameter. For this purpose, we use the data on 
J/ip production from the EMC collaboration |]T8|] in the limited z range 0.3 < z < 0.6. 



In this range, the colour- singlet contribution gives almost the entire cross-section |13 



We can therefore use these data to extract the effect of the smearing on the singlet 
wave-function. For (kx) = 0.7, 1.0 GeV, we find the values to be 0.98 GeV^ and 
1.41 GeV'^ respectively. The values of the non-perturbative octet parameters extracted 
for the three choices of (kx) are tabulated in Table 1. We see that the effect of the 
kx smearing on the parameters extracted from the large-p^ CDF data is very modest. 
Fig. 1 shows that the fits to the data when k^ smearing is included are very good and 
comparable in quality to the case (fcy) = 0. 



(kr) (GeV) 


Ml/ml + M2/3 (GeV)3 


M3 (GeV) 3 


0.0 
0.7 
1.0 


(3.14 ±0.58) X 10-2 
(2.82 ±0.47) X 10-2 
(2.35 ±0.39) X 10-2 


(1.26 ±0.33) X 10-2 
(1.35 ±0.30) X 10-2 
(1.50 ±0.29) X 10-2 



Table 1: The values of the colour-octet non-perturbative matrix elements determined 
from the CDF data on J/ ip production, for three choices of intrinsic initial state trans- 
verse momentum. 



Taking these fitted values of the parameters, we next consider inelastic J/ ip photo- 
production at HERA. For photoproduction, the formalism for the inclusion of smearing 
effects is similar to that used for the pp collisions described above, except that now 
we have only one hadron in the initial state. We take the same choice of parton dis- 
tributions, scales etc. as used in the Tevatron fits. We compute the z distribution 
for a photon-proton centre-of-mass energy y/s^p = 100 GeV, imposing a transverse 
momentum cut > 1 GeV. Again we present results for (kx) = 0, 0.7, 1.0 GeV. For 
each choice of (kx), the values of the octet non-perturbative matrix elements are taken 
from Table 1, and the singlet matrix elements are taken to be the same as that used 
in the Tevatron fits. 

With these inputs, we compute the z distribution for the three different values 
of {kx)- The results are compared with HERA data in Fig. 2. In the absence of 
smearing, (kr) = 0, we see that the colour-octet component makes a large contribution 
at z close to 1 which is not supported by the data. However the introduction of fc^ 
makes a substantial change to the octet contribution. Whereas the effect of /cT^-smearing 
is very small for large-pr production at the Tevatron, these effects are found to be very 
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Figure 2: The HERA data for da/dz (in nb) for J /ip production at J's^ = 100 GeV 



with pt > 1 GeV, compared with model predictions for three choices of the intrinsic 
transverse momentum distribution, namely (fcr) = 0, 0.7, 1.0 GeV. 



important for J/ip production at HERA. In particular, smearing signficantly reduces 
the size of the cross section (mainly due to the smearing function itself, but also due 
to the smaller octet matrix elements, see Table 1.) The z distribution also becomes 
flatter, in better agreement with the HERA data. In fact, excluding the highest two 
data points at 2; ~ 0.86 we find that the is minimised for (kx) ~ 0.7 GeV and the 
resulting description of the data is very good. It is perhaps not surprising that the 
highest-z data points are not exactly accounted for, since in this region other effects 
such as contamination from elastic scattering or breakdown of NRQCD factorisation 
may be important. 

In conclusion, while a direct comparison of the predictions of NRQCD with the 
z-dependence of the inelastic phot opro duct ion cross section for J/ip at HERA show a 
marked disagreement between the two, we argue that such a comparison is misleading. 
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The inelastic photoproduction process does not provide a clean test of NRQCD because 
in order to have a sizeable event rate, a very low p^-cut (~ 1 GeV) is necessary in 
the HERA experiments. At such low values of pt (and for values of z reasonably 
close to one), we find that the effect of fcr-smearing is important and that, indeed, for 
(fcr) ~ 0.7 GeV, the discrepancy between theory and experiment is no longer observed. 
On the other hand, the inclusion of smearing has a very modest effect on the large- 
ly J/i^ data from the Tevatron. Better tests of NRQCD may be obtained by studying 
other observables at the Tevatron itself, such as the study of the polarisation of the 
produced J/ip [|l9l or the production of other charmonium resonances [^0|, |2l| whose 
cross-section can be predicted in NRQCD. 
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